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Abstract 

A measurement of the W mass and width has been performed by the Delphi 
collaboration using the data collected during 1998. The data sample has an 
integrated luminosity of 155 pb _1 and an average centre-of-mass energy of 188.6 
GeV . Results are obtained by applying the method of direct reconstruction of 
the mass of the W from its decay products in both the W + W~ — > tvgofl and 
W + W~ — ► qq'qq' channels. The W mass result for the 1998 data set is 

M w = 80. 387 ± 0. 087 (stat) ± 0. 03A(sy st) ±0. 017 (LEP)±0.035(F SI) GeV/c 2 , 

where FSI represents the uncertainty due to final state interaction effects in 
the qq'qq' channel, and LEP represents that arising from the knowledge of the 
beam energy of the accelerator. Combining this result with those previously 
published by the Delphi collaboration gives the result 

M w = 80.359 ± 0.074(stat) ± 0.032(syst) ± 0M7(LEP) ± 0.033(FSI) GeV/c 2 . 
The combined value for the W width is 

T w = 2.266 ± 0.176(stat) ± 0.056(syst) ± 0.052(FSI) GeV/c 2 . 



(Accepted by Physics Letters B) 
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1 Introduction 

The W mass has been measured by the Delphi collaboration using the data collected 
during 1998. This data sample has allowed a significant improvement in the accuracy 
of the collaboration's W mass determination as the integrated luminosity is more than 
twice that on which previous Delphi results are based [0,0,0 • The W mass has also 
been determined by the other Lep collaborations [§] and at hadron colliders ||. Using 
the same reconstruction method as for the W mass, results on the direct measurement 
of the W width are also obtained in this paper, and can be compared with those of the 
other Lep collaborations and of the Cdf collaboration ||. 

Section || of this paper describes the characteristics of the 1998 data sample and of 
the event generators used in this analysis. 

The analysis is performed through the direct reconstruction of the mass of the W 
from its decay products in the W + W~ — > qq'qq' (fully-hadronic) and W + W~ — > iP^qq' 
(semi-leptonic) decay channels. The applied methods are described in Sec. [|, and have 
been refined from those in previous publications. Results are now reported for rV T qq' 
events, and in both the fully-hadronic and semi-leptonic channels improvements in the 
handling of initial-state radiation (ISR) have been included. 

The systematic error evaluation described in Sec. |] has increased in sophistication 
from that previously reported. New techniques have been applied, such as the use of 
mixed Lorentz boosted Z's, and a wider range of higher precision simulation studies have 
been performed. 

The results of this analysis are reported in Sec. |5|, and are combined with the previous 
Delphi results. 

2 Data and Simulation Samples 
2.1 Data 

A detailed description of the Delphi apparatus and its performance can be found 
in [0] []. In 1998 the detector was used to record data at the Z peak and at a nominal 
centre-of-mass energy of 189 GeV. 

The Z peak data were recorded before (1.8 pb _1 ) and towards the end (0.6 pb _1 ) of 
the high energy data, thus facilitating checks of detector stability. These data were the 
principal sample used for calibration and alignment of the detector and, in this analysis, 
assist the study of systematic uncertainties. 

The average centre-of-mass collision energy for the high energy data was 188.6 GeV. 
The luminosity weighted r.m.s. of this value, assessed on a fill by fill basis, was 50 MeV. In 
the data sample considered for analysis all the detectors essential for this measurement 
were required to be fully efficient; the operation of the central tracking detectors was 
important for all decay channels, in the iveqq' analysis stricter requirements than in the 
qq'qq' channel were placed on the electromagnetic calorimeters. The selected samples 
correspond to integrated luminosities of 152.9 pb _1 for the iPeqq' analysis and 157.4 pb _1 
for the qq'qq' decay channel. 



1 The co-ordinate system used has the z-axis parallel to the electron beam, and the polar angle calculated with respect 
to this axis. 
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2.2 Simulation 

The response of the detector to various physical processes was modelled using the sim- 
ulation program DELSIM @, which includes modelling of the resolution, granularity and 
efficiency of the detector components. In addition, detector correction factors, described 
in Sec. f|, were included to improve the description of jets, electrons and muons. For 
systematic uncertainty studies a fast simulation program, relying on a relatively simple 
set of smearing and efficiency parametrisations, was also used. 

WW events and all other four-fermion processes were produced using the event gen- 
erator EXCALIBUR ||, with initial-state radiation described using the QEDPS program ||. 
The W mass (M w ) and width (T w ) definition used throughout this paper correspond to 
a W propagator with an s-dependent width. The background process e + e~ — > qqi'j) 
was simulated with the PYTHIA 5.7 [TO] event generator. Two photon backgrounds were 



studied using the TWOGAM generator ||11|| . The fragmentation of all events was performed 
using JETSET 7.4 [Kj tuned to the DELPHI LepI data [|1^]. Systematic error checks were 
performed using other generators and variations in the fragmentation tuning as described 
in Sec. |j. The systematics section also reports Delphi results for the common Lep event 
samples produced in the context of the Lep WW Workshops: the event generation was 



performed by ALEPH using KORALW 1.21 [13], and these events were then passed through 
the DELSIM program. 



3 Analysis Method 

3.1 Semi-Leptonic Decay Channel 

The fitting procedure presented here is a development of that described in || for the 
e!7 e qq' and //z7 M qq' channels, where the fitting function now includes a description of the 
ISR spectrum of WW events. In addition, we present here an analysis of the rV T qq' 
channel and the event selection in all semi-leptonic channels is now based on a neural 
network. 



3.1.1 Event Selection 
Lepton Identification 

Charged particles were identified as muons if they were associated with a hit in the 
muon chambers, or had an energy deposit in the hadron calorimeter that was consistent 
with a minimum ionising particle. Muon identification was performed in the polar angle 
range between 10° and 170°. 

Electron identification was performed in the polar angle range between 12° and 168° by 
selecting charged particles with a characteristic energy deposition in the electromagnetic 
calorimeters. In the central region of the detector covered by the HPC electromagnetic 
calorimeter, electrons were selected using the energy over momentum (E/p) ratio of the 
candidate. For lower energy candidates (below 30 GeV) this was supplemented by selec- 
tion criteria on the shape of the calorimetric shower and a more stringent comparison of 
the track extrapolation and calorimetric shower positions, while for electrons of higher 
energy, negligible energy deposition in the hadron calorimeter was required. In the po- 
lar angle region below 36° and above 144°, where the momentum resolution is poorer, 
tracks associated to electromagnetic energy showers above 8 GeV and negligible hadron 
calorimeter energy deposits were considered as electrons. 
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Tau candidates were obtained by clustering the events into a three-jet configuration 



using the LUCLUS [14] algorithm. The jet with the lowest charged multiplicity was chosen 
as the tau candidate. As the tau lepton predominantly decays into a final state with one 
or three charged particles, only jets containing between one and four charged tracks were 
selected. 

Selection 



The event selection was based upon a multi-layer perceptron neural network [15||. The 



network was optimised separately for ez7 e qq', /il^qq', rz7 T qq' candidates containing only 
one charged particle and rV T qq' candidates with several charged particles. 

Having removed the lepton candidate in ez7 e qq' and /xX^qq' events, the LUCLUS jet 
clusterization algorithm (with a djdn of 7.5 GeV/c) was used to cluster the remaining 
particles. Events containing more than three jets were re-clustered, forcing them into a 
three-jet configuration. rz7 T qq' events were clustered as the tau candidate and a two-jet 
system. The events were reconstructed using a constrained fit imposing conservation of 
four- momentum and equality of the two W masses in the event. As the energy of the tau 
lepton is unknown, due to the emission of at least one neutrino in its decay, the mass in 
the rlAj-qq' channel is entirely determined by the jet system. 

The neural network relied upon the characteristic event properties in each decay chan- 
nel. The input variables included lepton momentum, polar angle of the missing momen- 
tum, and the isolation of the lepton candidate from the hadronic system of the event. The 
electron and muon identification was obtained as a strong or loose tag and the network 
was optimised separately in each decay channel for these two categories of identification. 

The network was tuned on samples of signal and background simulation events, and 
its performance estimated on independent samples of events. After applying a cut on the 
network output the selection performance is as shown in Tab. |l|. Events that passed the 
cut in the muon channel were selected, the remaining events were considered as electron 
channel candidates and, if they were again rejected, were then analysed under the tau 
channel hypothesis. The neural network output in the r77 r qq' channel is shown in Fig. [I]. 

The tau selection sample contains a significant proportion of other semi-leptonic de- 
cays: the composition was estimated from simulation to be 67% T77 T qq', 19% ez7 e qq', 5% 
/iz7 M qq' with the remaining fraction dominated by qqi'j) background events. This corre- 
sponds to a 39% selection efficiency for rz7 T qq' events. Further information on a selection 



of W + W~ events in Delphi with a similar performance is available in \W . 

The fraction of semi-leptonic WW events in the sample was extracted from simulation 
as a function of the neural network output: this is referred to below as the event purity 
P e . This feature is particularly useful for the tau selection, where the proportion of 
background events is highest. 



3.1.2 Likelihood Function 

The following likelihood function was evaluated for each selected events with a recon- 
structed mass in the range 68 — 92 GeV/c 2 : 

4(M W , T w ) = P e ■ S"(m-f u , a*, M w , T w ) + (1 - P e ) ■ B(m-f u ), (1) 

where P e is the event purity, discussed above, S" is the signal function that describes the 
reconstructed mass distribution of the semi-leptonic W decays, and B is used to describe 
background processes. The reconstructed event mass vrJ %t and its estimated error o^ 1 
are both obtained from the constrained fit. The distribution of background events is 
extracted from simulation SIS db function of m* %t . 
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The signal function S" is defined in terms of S,S' as discussed below. The function S 
relies on the convolution of three components, using x and m as the dummy integration 
variables: 

S{m ftt ,a flt ,M w ,T w ) = / dm G[m flt -m,a flt } / dx BW PS [m-(\-x), M w R ISR (x)]. 

Jo Jo 

(2) 

BWps is a phase-space corrected relativistic Breit-Wigner distribution (representing 
the W mass distribution) which is convoluted with the Gaussian function G describing 
the detector resolution. The width of the Gaussian depends upon the reconstructed mass 
error obtained in the constrained fit for that event. Details of the BWps and G terms 
are given in 0. A recent addition to the analysis is the description of the ISR spectrum, 
which is parametrised as 

RiSRixJ = Px^, (3) 

where x 7 is the ratio of the photon energy to the centre-of-mass energy and /3 is calculated 
from the electromagnetic constant (a), the centre-of-mass energy squared (s) and the 
electron mass (m e ): 

In 

P = — [\og(s/m 2 e )-l]. (4) 

7T 

Including this ISR term decreases the bias on the fitted W mass by approximately 
400 MeV and improves the expected error by 2 ± 1%. 

The event selection contains a significant fraction of r77 T qq' events in the electron and 
muon channel samples, and of e77 e qq' events in the tau sample (see Tab. [I]). In the tau 
channel the mass of the event is determined from the jet system. The behaviour of true 
TV T qq' and ez7 e qq' events in this fit are found to be similar. However, in the electron and 
muon channel samples the behaviour of the rV T qc[ events is somewhat different to that 
of the ez7 e qq', fiV^qq' events. The TV T qq' events have a worse mass resolution and a small 
negative bias on the mass. The fraction of tau events, which have been wrongly classified 
and are contained in the electron and muon channel samples, has been parametrised in 
bins of the lepton energy and the measured missing mass. This event impurity Pr e was 
then taken into account in the likelihood function for the electron and muon samples, by 
defining the signal function S" as 

S" = (l-Pr e )-S + Pr e -S', (5) 

where S' is analogous to S, but with the width of the Gaussian resolution function 
increased according to simulation studies. All remaining biases in the analysis due to 
using this approximate likelihood description are corrected for in the calibration procedure 
as described in Sec. ^3| . 

3.2 Fully-Hadronic Decay Channel 

The analysis of the fully-hadronic channel was based on that applied in 0|. However, 
the implementation now relies on kinematic fits with four rather than six constraints and 
includes a new ISR treatment. 
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3.2.1 Event Selection 

A sample of hadronic events was selected by requiring more than 13 charged particles 
and a total visible energy exceeding 1.15 -E BEA m- 

qq(j) events were suppressed by demanding an effective centre-of-mass energy ]T7| , 
after ISR emission, of greater than 161 GeV. The algorithm for assessing the e + e~ 
collision energy considers both the emission of unobserved ISR photons in the beam-pipe 
and photon candidates detected in the electromagnetic calorimeters. 



The DURHAM jet clustering algorithm [18] with y cut of 0.002 was applied to the event. 



If any of the resulting jets contained less than three particles or had an invariant mass 
smaller than 1 GeV/c 2 , clustering was continued to a higher value of y cu t. Events with 
less than four jets were then rejected, while events containing six or more jets were re- 
clustered into five objects representing four quarks plus one hard gluon jet. 

Events containing b-quarks were rejected using the Delphi b-tag package [Tj|, this 



removes 17% of ZZ events and 6% of qq(j) background while reducing the signal efficiency 
by only 0.2% . A four constrained kinematic fit |5] was applied to the remaining events, 
enforcing conservation of energy and momentum. 

A variable to discriminate between qq events with hard gluon radiation and signal 
events was constructed. This compound variable relied upon the fitted jet energies and 
the inter-jet angles. The expected fraction of qq'qq' events in the selected sample, the 
event purity, was parametrised as a function of this variable. Events with an estimated 
purity below 25% were rejected. 

The performance of the event selection is shown in Tab. Further information on a 



selection of W + W events in Delphi with a similar performance is available in |16| . 



3.2.2 Likelihood Function 

For each of the selected events an event likelihood was constructed : 

£ e (M w , T w ) = JjYl w ^ ' PiA m xi m v) • (6) 

i 

Pf f - S , (m x ,?7i w ,M w ,r w ) + (1 - P e J s ) ■ B] dm x dm y . 

As in |3| the signal function ^(m^mj^MwTw) consists of Breit-Wigner terms for the 
WW and the ZZ contribution and a phase space correction factor. A flat distribution 
B accounts for background processes and wrong jet pairings in the signal events. Both 
S and B were normalised to 1 over the integration area. The fraction of the signal and 
background likelihoods used for each event depend upon the event purity P^^ ' . This 
purity was parametrised as a function of a discriminating variable as described above. 

The sum J2i w i,e • Vi,e{. m xi m y) is a weighted sum of the probability densities p i;f , that 
the event e corresponds to two heavy objects with mass m x and m y . 

The probability density Pi je (m x , m v ) oc exp — Jjn xi m v ) was determined for all jet 
pairings (three possibilities for a four-jet event and ten for a five-jet event) and with three 
different clustering algorithms (DURHAM @, CAM JET [gOf and DICLUS f2l|). The relative 



probabilities w ifi that the corresponding jet pairing was the correct one were estimated us- 
ing jet charge information and the transverse momentum of the gluon candidate (see ||). 
The three jet clustering algorithms were given the same weight. 

The probability was calculated using a Gaussian approximation for the x 2 '- 



Xle(m x , my) « xlc + (m - m Rt ) T \-\m - m fit ) (7) 
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with 

V 



< WpP- ) , m = ( M and m fit = { < 



The masses m~*, m* 4 , their errors cr ma . and a my and the correlation between them, p xy , 
are taken from a four constrained kinematic fit. When the x\c * s larger than the number 
of degrees of freedom (NDF = 4), the xl e ( m x, m y) is rescaled with a factor NDF/xlc 
in order to compensate for non-Gaussian resolution effects. This procedure decreases the 
computing time taken by an order of magnitude compared with the full six constrained 
fit ||, while resulting in only a minimal reduction in the W mass precision obtained 
(2±1%). 

A new feature of this analysis is a treatment of events under the collinear ISR hypoth- 
esis. A kinematic fit was performed with modified constraints to simulate the emission 
of an ISR photon of momentum p z inside the beam pipe: 

^objects 

£ {E, Px ,p y , Pz ) i = {^rs-\ V %^pf). (8) 

i=l 

The probability that the missing momentum in the z direction was indeed due to an 
unseen ISR photon was extracted from the simulation as a function of |pf t |/o" Pz , where 
a Pz is the estimated error on the fitted z momentum component; only events with this 
ratio greater than 1.5 were treated with the mechanism described below. 

Additional probability density Pi >e terms were then included in the likelihood sum for 
these events, with a relative weight factor derived from the ISR hypothesis probability. An 
example of the effect of including the ISR hypothesis is shown in Fig. |2|. This treatment 
was applied to 16% of the events and resulted in an improvement of the expected W mass 
error for these events of 15%. 



3.3 Mass and Width Extraction 

The distribution of the reconstructed invariant masses of the selected events after 
applying a kinematic fit, imposing four- momentum conservation and the equality of the 
two di-jet masses, are shown in Fig. |3|. This plot is provided for illustrative purposes 
only, the mass and width fitting procedure is described below. 

The combined likelihood of the data can be obtained from the product of the event 
likelihoods described above. The W mass and width were extracted from maximum 
likelihood fits. The W mass fit is performed assuming the standard model value for the 
W width. The W width was obtained assuming a mass of 80.35 GeV/c 2 . The correlation 
between M w and r w was found to have a negligible impact on the extracted width value. 

The mass and width analyses have been calibrated separately in each of the decay 
channels (qq'qq^eI7 e qq^/iI^ M qq^rI7 r qq'). The biases of the analyses were estimated by 
re-weighting generated simulation samples to obtain the fitted mass and width values. 
The re-weighting was performed using the extracted matrix element of the EXCALIBUR 
generator. The linearity of the mass analysis was estimated using independent simulation 
samples generated at three W mass values, while the re-weighting procedure was used 
for the width analysis. The analyses were corrected with the calibration results, and the 
statistical error on the bias is included in the systematic error. 

The analyses were checked by performing fits to a large number of samples of simulation 
events. Each sample was comprised of a mixture of signal and background simulation 
events to represent the expected distribution in the data. The pull distribution ms%t a ™ w 
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was demonstrated to be compatible with a Gaussian of width one to an accuracy of better 
than 1%. The mean expected statistical error in the W mass was 262 MeV/c 2 for ez7 e qq', 
203 MeV/c 2 for /il/^qq', 311 MeV/c 2 for rz7 T qq' and 104 MeV/c 2 for the qq'qq' channel. 



4 Systematic Uncertainties 

The sources of systematic error that have been considered for the W mass and width 
determinations are described in the subsections below. The results of these studies are 
summarised in Tabs. || and |[ 



4.1 Calibration 

The accuracy with which the bias of the analysis can be determined is limited by 
the size of the generated simulation samples. Sufficient events were generated to limit 
this error to less than 10% of the statistical error in any given channel. The calibration 
procedure is described in Sec. |3~5. 



4.2 Detector Effects 

The data taken at the Z peak were used to study, and limit, possible errors in the 
detector simulation model. 

Muon studies were performed on a selected sample of Z — > events. From the di- 

muon sample corrections to the inverse momentum scale, 1/p, were calculated separately 
for positive and negative muons as a function of the lepton polar angle. The systematic 
error on this correction was estimated by varying it by half of its value. The momentum 
resolution (typically 0.001 in 1/p) was found to be slightly better in simulation than in 
the data (a maximum difference of 10%). This was corrected by smearing the simulation 
with a Gaussian. An extra smearing of 0.0005 in 1/p was used to estimate the systematic 
error coming from this correction. The combined systematic error from these corrections 
is quoted for the /w^qq' channel as the lepton correction systematic error in Tabs. [| and 

s 

The correction of the energy scale of electrons was determined from Bhabha events at 
the Z peak in different polar angle regions. The residual systematic error on this absolute 
energy scale was estimated to be 0.5%. In each of these polar angle regions, the energy 
resolution of simulation events was degraded by applying a Gaussian smearing, and the 
residual error on this smearing was estimated to be 1%. The dependence of the energy 
calibration as a function of the electron energy was checked using low energy electrons 
from Compton events at the Z peak, and high energy electrons from radiative Bhabha 
scattering at high centre-of-mass energy. In these cases the true energy of the lepton 
was deduced from 3-body kinematics using only the angular information and assuming 
that the unseen particle was along the beam axis. The absolute energy calibration was 
found to be compatible with requiring no additional corrections in all energy ranges. The 
systematic error coming from this source was estimated assuming a 1% change of slope 
in the energy calibration in the range of interest (between 25 and 70 GeV). The lepton 
correction uncertainty for the ez7 e qq' channel, quoted in Tabs. || and |9], is the quadratic 
sum of the errors from these three error sources. 

The lepton correction error is not quoted for the r77 T qq' channel as the reconstructed 
tau lepton carries no information for the reconstructed W mass. 
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Jet energies were studied in Z — > qq events as a function of the polar angle and 
reconstructed energy of the jet. The comparison of data and simulation showed agreement 
within a band of ±2% over most of the Delphi detector's angular coverage and an overall 
uncertainty of 1% was estimated. The simulation was smeared by a Gaussian function to 
improve the description of the observed energy spread in the data and the residual error 
on this smearing was estimated to be ±4%. The dependence of the energy calibration 
as a function of the jet energy was checked using low energy jets from qq + gluon events 
at the Z peak and high energy jets from radiative Z's at high centre-of-mass energy. 
The true jet energy was estimated from 3-body kinematics; in the radiative events the 
unseen photon was assumed to be along the beam axis. No additional energy calibration 
slope was necessary over the relevant energy range (25 to 75 GeV) and a 1% change in 
slope was used to calculate the systematic error. A study of the acollinearity of jets in 
Z — > qq events was performed and an appropriate smearing to the simulation of the jet 
angular direction was estimated. A systematic error was estimated by applying an extra 
5 mrad angular smearing. The jet correction uncertainty, quoted in Tabs. |8] and ^, is the 
quadratic sum of these four errors. 

A possible source of angular distortion in Delphi is the uncertainty on the length 
to width ratio of the detector. The detector is aligned relative to the vertex detector, 
the largest uncertainty being the radius of this detector which is known to a precision 
of ±0.1%. This error is listed in the Mw systematics table under aspect ratio. The 
corresponding uncertainty for the width measurement is less than 10 MeV. 



4.3 Background Description 

The background level was changed by ±10% in the simulation; this easily covers the 
expected uncertainty in the accepted cross-section, as discussed in |fT6|| . The dominant 
background source, Z — > qql'j), was generated using both JETSET and HERWIG fragmen- 
tation models and mass fits performed using both these samples. The 4-jet rate in data 
and simulation was also studied using events collected at the Z peak. It is concluded that 
the background description is a relatively small component of the systematic error for 
the mass measurement, as shown in Tab. [8] and somewhat more important for the width 
(see Tab. |). 



4.4 Fragmentation 

A study of the possible effects on M w and T w of the simulation of the event fragmen- 
tation was performed, the results are provided in Tabs. and |3| . 

Events were produced using a fast simulation package where the value of Aq C d and 
<j q were changed with respect to the standard Delphi JETSET tuned values |12| by twice 
their estimated errors Q. The estimated errors are ±0.018 GeV for Aqcd and ±0.007 GeV 
for a q . 

The event samples generated by Aleph in the context of the LEP WW workshop are 
used to compare results from the HERWIG and JETSET fragmentation models. The HERWIG 



events were produced with a recent tuning |22] which provides a better description of the 
data than previous HERWIG versions. 

The results on the Delphi tuning of JETSET are also compared with those of the 
Aleph collaboration. This comparison cross-checks several effects as the Aleph events 

2 The dominant systematic error components of the tuning uncertainties were estimated from a comparison of fits with 
a range of input data distributions 
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were produced with a different generator (KORALW) from the Delphi events which includes 
a different ISR and final-state radiation (FSR) treatment. 

As the results in Tabs. ^| and |3| are all compatible with zero, we quote a systematic 
error from fragmentation reflecting twice the statistical precision of the JETSET tuning 
parameter studies in the combined semi-leptonic and fully-hadronic channels. 

4.5 Mixed Lorentz Boosted Z's 

The agreement between data and simulation can be cross-checked using the method 
of Mixed Lorentz Boosted Z's. Z events were selected from the Z peak data sample 
collected during 1998 and the corresponding simulation sample. Through a suitable 
choice of Lorentz boost and superimposing two Z events a WW event may be emulated. 
The angular distribution of the Z events used was chosen to match that expected in WW 
events. 

This technique has been applied to the qq'qq' mass measurement and the hadronic 
system in semi-leptonic events. The differences between data and simulation obtained are 
2 ± 2 MeV/c 2 and 3 ± 10 MeV/c 2 in the two topologies respectively, where the errors are 
statistical. These results reinforce the view that the quoted systematic errors for detector 
and fragmentation effects are conservative. A study of possible intrinsic uncertainties of 
the MLBZ method j23| in the qq'qq' channel estimates the accuracy of this technique to 
be 5 MeV/c 2 , and demonstrates excellent agreement between data and simulation as a 
function of relevant event variables. 



4.6 ISR 

A comparison of two independent models of ISR was performed. WW events produced 
with the EXCALIBUR generator were re-weighted as a function of the total ISR energy in 
the event. The weights were obtained for each WW decay channel from a generator level 
study of the KORALW [T3] ISR treatment (based on the YFS exponentiation approach) 
and that of the standard algorithm used in Delphi QEDPS (based on a parton shower 
approach). The results of the W mass and width fits are shown in Tabs. f| and |5|. We 
conservatively choose to quote the largest deviation observed in any of the channels. 



4.7 LEP Beam Energy 

The average Lep beam energy at Delphi is evaluated by the energy working group 



at 15 minute intervals of running or after significant changes in the beam energy. 



The measured centre-of-mass energy is imposed as a constraint in the kinematic fit, and 
hence the relative error on the beam energy translates to the same fractional error on the 
W mass determination. The spread of energies of the electrons and positrons in the LEP 
beams was found to have a negligible impact on the mass and width measurements. 



4.8 Bose-Einstein Correlations 

Bose-Einstein statistics dictate that the production amplitude for final state particles 
should be symmetrical under the exchange of identical bosons. The omission of these 
correlations between particles from different W bosons in our standard simulation could 
lead to a systematic error on the W mass and width measurement in the qq'qq' channel. 
A clear picture has yet to emerge from the experimental study of this phenomenon 
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To evaluate the possible size of the effect on the mass and width measurements we 
have considered several phenomenological models; the results are given in Tab. |]. The 
relevant value for the systematic uncertainty is the difference between the shifts obtained 
from Bose-Einstein correlations inside individual W's and that between W's. The models 
used are: 



The LUB0EI algorithm [25| in JETSET changes the momentum vectors of identical 
final state bosons to model the two particle correlation function, and then offers a 
range of options to restore energy and momentum conservation. 
Global re-weighting aims to reproduce the enhancement of identical bosons close 
in phase-space by giving weights to events. This procedure does ensure energy and 
momentum conservation, but may adversely affect other event distributions. DELPHI 
results were obtained using the Kartvelishvili/Kvatadze/M0ller (KKM) re-weighting 



scheme |27j and reported on in ||. 

In H we also reported results for a study based on a modification of the JETSET 
fragmentation model introducing quantum mechanical 2-particle and 3-particle in- 



terference (ST) for identical bosons using a local re-weighting technique p8| 



The version of LUB0EI model studied here is BE 32 , in which a local energy and mo- 
mentum conservation procedure is applied Six sets of fully simulated Z events were 



generated at a range of values of correlation strength (A) and radius (r), and the four 
momentum difference, Q, between all selected same charge particle pairs was calculated. 
An interpolation was performed on the basis of the Q distribution to obtain A = 1.35 
and r = 0.6 fm (PARJ(93) = 0.34 GeV) which provide the optimal description of the 
Delphi Z peak data. The simulated event samples were produced with the standard 
DELPHI JETSET tuning. The W mass and width shifts were evaluated using samples 
produced with fast simulation. 

4.9 Colour Reconnection 

The hadronization of two W bosons in the qq'qq' channel may not occur independently. 
The colour flow between the W's may lead to a shift in the measured W mass. Although 



experimental work is progressing | 29fl , a suitable sensitivity has not yet been reached by 
the measurements to limit the effect on Mw . 

In the previous publication we reported results with the ARIADNE colour reconnec- 
tion model, the results are repeated in Tab. [F[ However, we note that the current version 
of the ARIADNE 2 model is disfavoured by LepI data and that we do not consider 
the ARIADNE 3 model for the systematic error assessment as it allows perturbative phase 



reconnection where calculations have shown the effect to be small 



The Sjostrand/Khoze models of colour reconnection are available in the JETSET frame- 
work. We have used the SKI and SK2 models. Results for the SKI model are quoted for 
30% of reconnected events, this is the same reconnected fraction as in the SK2 model. 
Events were produced using DELSIM and a fast detector simulation. By processing the 
same event sample through the full detector simulation and the fast simulation, the 
reliability of the fast detector simulation for this study was clearly demonstrated. The 
simulated event samples were produced with the standard Delphi JETSET tuning. 

However, we report that the standard implementation of SKI shows numerical in- 
stabilities that reduce the reliability of the model used. In this model the reconnection 
probability is a function of the string overlap. This overlap is calculated by numerical 
integration through sampling. The accuracy of this calculation is improved by using 1000 
sampling points rather than the default value of 100. In addition, we also report results 
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from a more efficient sampling of the string overlap, in which the sampling is performed 
along strings taking into account their life-time and the total overlap is calculated as the 
sum of the overlap of pairs of string pieces. 

The observed W mass and width shifts are given in Tab. |^ and in Fig. [| the observed 
W mass and width shifts are shown as functions of the reconnected fraction of events in 
the SKI model. 

We conservatively choose to quote as a systematic error the largest effect observed in 
our studies; this is approximately 50 MeV for both the mass and width analyses. 



4.10 Correlations 

The components of the systematic error arising from the jet energy scale corrections, 
aspect ratio, ISR, fragmentation and LEP beam energy are taken as correlated between 
the analyses in the different decay channels for the 189 GeV data. The background 
description and lepton modelling uncertainties are treated as uncorrelated between WW 
decay channels. 

In Sec. f).2\ a combination with the previously published Delphi results is performed. 



The LEP energy correlation matrix is used [^] in this process. The calibration statistics 
error is uncorrelated between years, while all other systematic errors are conservatively 
assumed to be fully correlated between years. 



5 Results 



5.1 189 GeV Results 

The W mass and width results of the analyses described in this paper are presented in 
Tabs. |10]and[ll]. The error is divided into its statistical component, indicated (stat), the 
main systematic component (syst) and the systematic from the beam energy uncertainty 
(LEP). In the qq'qq' channel an error from final state interaction effects (FSI) is also 
included. 

The semi-leptonic and fully-hadronic W mass results are combined and the following 
result obtained: 

M w = 80.387 ± 0.087(stat) ± 0.034(syst) ± 0M7(LEP) ± 0.035(FSJ) GeV/c 2 . 
The combined result for the W width is : 

T w = 2.205 ± 0.195(stat) ± 0.059(syst) ± 0M7(FSI) GeV/c 2 . 



5.2 Combined DELPHI Results 

These results are combined with the previously published Delphi results on the W 
mass determination by direct reconstruction 0j3|] and the result from the threshold cross- 
section measurement [|TJ to obtain: 

M w = 80.359 ± 0.074(stat) ± QM2(syst) ± 0M7(LEP) ± 0.033(FSJ) GeV/c 2 . 

This combination has a x 2 of 4.5 with 6 degrees of freedom. The combined result on the 
W mass from the qq'qq' channel alone is : 

M w = 80.369 ± 0.091(stat) ± 0.029(syst) ± 0M7(LEP) ± 0.056(FSJ) GeV/c 2 
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and for the £z/^qq' channel alone is : 

M w = 80.327 ± 0.128(stat) ± 0M5(syst) ± 0M7(LEP) GeV/c 2 . 

The mass difference between the measurements for the fully-hadronic and semi-leptonic 
channels, AMw(qq'qq' — ^qq'), nas been determined: 

AM w (qq'qq' - fl7/qq') = 39 ± 159 MeV/c 2 . 

A significant non-zero value for AMw could indicate that FSI effects are biasing the 
value of Mw determined from the qq'qq' events. As AMw is primarily of interest as a 
check of the possible effects of final state interactions, the errors from FSI effects are set 
to zero in this determination: all other errors and correlations were as described above. 

The combination for the W width with the 1997 data [Q yields : 

T w = 2.266 ± 0.176(stat) ± 0.056(syst) ± 0.052(FSJ) GeV/c 2 . 

Acknowledgements 

We are greatly indebted to our technical collaborators, to the members of the CERN- 
SL Division for the excellent performance of the LEP collider, the LEP energy working 
group for the beam energy estimation and to the funding agencies for their support in 
building and operating the DELPHI detector. 
We acknowledge in particular the support of 

Austrian Federal Ministry of Science and Traffics, GZ 616.364/2-III/2a/98, 
FNRS-FWO, Flanders Institute to encourage scientific and technological research in the 
industry (IWT), Belgium, 

FINEP, CNPq, CAPES, FUJB and FAPERJ, Brazil, 

Czech Ministry of Industry and Trade, GA CR 202/96/0450 and GA AVCR A10 10521, 

Danish Natural Research Council, 

Commission of the European Communities (DG XII), 

Direction des Sciences de la Matiere, CEA, France, 

Bundesministerium fur Bildung, Wissenschaft, Forschung und Techno logie, Germany, 
General Secretariat for Research and Technology, Greece, 

National Science Foundation (NWO) and Foundation for Research on Matter (FOM), 
The Netherlands, 
Norwegian Research Council, 

State Committee for Scientific Research, Poland, 2P03B06015, 2P03B11116 and 
SPUB/P03/DZ3/99, 

JNICT- Junta Nacional de Investigagao Cientifica e Tecnologica, Portugal, 
Vedecka grantova agentura MS SR, Slovakia, Nr. 95/5195/134, 
Ministry of Science and Technology of the Republic of Slovenia, 
CICYT, Spain, AEN96-1661 and AEN96-1681, 
The Swedish Natural Science Research Council, 
Particle Physics and Astronomy Research Council, UK, 
Department of Energy, USA, DE-FG02-94ER40817. 



13 



References 

[1] DELPHI Collaboration, P. Abreu et al, Phys. Lett. B397(1997) 158 
[2] DELPHI Collaboration, P. Abreu et al, Eur. Phys. J. C2(1998) 581 
[3] DELPHI Collaboration, P. Abreu et al, Phys. Lett. B462(1999) 410 
[4] ALEPH Collaboration, R. Barate et al, Eur. Phys. J. CC17(2000) 241 
L3 Collaboration, M. Acciarri et al, Phys. Lett. B454(1999) 386 
OPAL Collaboration, G. Abbiendi et al, hep-ex/0009018 (2000) 
[5] CDF Collaboration, T. Affolder et al, hep-ex/0007044 (2000), to be published in 
Phys. Rev. D 

DO Collaboration, B. Abbott et al, Phys. Rev. Lett. 80(1998) 3008 

DO Collaboration, B. Abbott et al, Phys. Rev. Lett. 84(2000) 222 
[6] CDF Collaboration, T. Affolder et al, Phys. Rev. Lett. 85(2000) 3347 
[7] DELPHI Collaboration, P. Aarnio et al, Nucl. Instr. and Meth. A303(1991) 233 

DELPHI Collaboration, P. Abreu et al, Nucl. Instr. and Meth. A378(1996) 57 
[8] F.A. Berends, R. Pittau, R. Kleiss, Comput. Phys. Commun. 85(1995) 437 
[9] Y. Kurihara, J. Fujimoto, T. Munehisha, Y. Shimizu, Progress of Theoretical Physics 

Vol. 96 (1996) 1223 
[10] T. Sjostrand, Comput. Phys. Commun. 82(1994) 74 

[11] S. Nova, A. Olchevski and T. Todorov in Physics at LEP2, CERN 96-01 Vol.2 (1996) 
224 

[12] DELPHI Collaboration, P. Abreu et al, Zeit. Phys. C73(1996) 11 
[13] S. Jadach, W. Placzek, M. Skrzypek,B.F.L. Ward, Z. Was, Comput. Phys. Commun. 
119(1999) 272 

[14] T. Sjostrand, PYTHIA 5.7 and JETSET 7.4: Physics and manual, CERN-TH-7112- 
93-REV (1995) 

[15] Code kindly provided by J. Schwindling and B. Mansoulie 

[16] DELPHI Collaboration, P. Abreu et al, Phys. Lett. B479(2000) 89 

[17] P. Abreu et al., Nucl. Instr. and Meth. A427(1999) 487 

[18] S. Catani, Yu.L. Dokshitzer, M. Olsson, G. Turnock, B.R. Webber, Phys. Lett. 
B269(1991) 432 

N. Brown, W. Stirling, Zeit. Phys. C53(1992) 629 
[19] G. Borisov, Nucl. Instr. and Meth. A417(1998) 384 

DELPHI Collaboration, P. Abreu et al, Eur. Phys. J. C10(1999) 415 
[20] Yu.L. Dokshitzer, CD. Leder, S. Moretti, B.R. Webber, JHEP 08 (1997) 001 
[21] L. LSnnblad, Zeit. Phys. C58(1993) 471 
[22] G.Rudolph, ALEPH HERWIG tuning - private communication 
[23] N.Kjaer, M.Mulders, CERN-OPEN-2001-026 (2001) 
[24] LEP Energy Working Group, Eur. Phys. J. Cll(1999) 573 

LEP Energy Working Group, LEP Energy working group note 99/01 (1999) 
[25] A. Valassi, Bose-Einstein correlations in W decays, hep-ex/0009039 (2000), to be 

published in the proceedings of the 30th International Conference On High-Energy 

Physics, 27 Jul - 2 Aug 2000, Osaka, Japan 
[26] L. LSnnblad, T. SjSstrand Eur. Phys. J. C2(1998) 165 
[27] V. Kartvelishvili, R. Kvatadze, R. M0ller, Phys. Lett. B408(1997) 331 
[28] S. Todorova, J. Rames, IReS-97-29 PRA-HEP-97-16, hep-ph/9710280 (1997) 
[29] P. de Jong, Color reconnection in W decays, hep-ex/0103018 (2001), to be published 

in the proceedings of the 30th International Conference On High-Energy Physics, 27 

Jul - 2 Aug 2000, Osaka, Japan 



14 



[30] OPAL Collaboration, G. Abbiendi et ai, Eur. Phys. J. Cll(1999) 217 
[31] T. Sjostrand, V. Khoze, Zeit. Phys. C62(1994) 281 
T. Sjostrand, V. Khoze, Phys. Rev. Lett. 72(1994) 28 



15 







Event Selection 




Event Type 


ez/ c qq' 


/iX^qq' 


ri/ T qq' 




qq'qq' 

11 11 




259.5 


0.3 


43.6 


303.4 


4.2 




0.5 


319.6 


10.8 


330.9 


2.5 


rz7 r qq' 


12.1 


13.4 


155.4 


180.9 


8.0 


qq'qq' 


0.2 


0.2 


1.7 


2.1 


1112.7 


Other 4f 


3.4 


3.0 


4.4 


10.8 


0.0 


55(7) and other 2f 


8.0 


0.7 


15.8 


24.5 


346.9 


Total 


283.7 


337.2 


231.7 


852.6 


1474.2 


Data 


244 


307 


236 


787 


1481 



Table 1: Number of selected events from the 189 GeV data sample, and the corresponding 
number of expected events from the simulation. Column four is the sum of the three 
previous columns. 



M w Fragmentation study (MeV/c 2 ) 



Event Type 


ez7 e qq' 


/iz7 M qq' 


rz7 r qq' 


qq'qq' 


Aqcd - 2a 
Aqcd + 2(T 

(T q — 2<J 

(Jq + 2a 


0± 14 
-4 ±14 
0± 14 
0± 14 


+9± 16 
±10 ± 16 
±30 ± 16 
±30 ± 16 


±30 ± 30 
-10 ±30 
-3 ±30 
±27 ±30 


±3 ±6 
-7 ±6 
-9 ±6 
-3 ±6 


Aleph herwig - Aleph jetset 
Aleph jetset 


75 ±52 
15 ±51 


8 ±42 
-17 ±42 


106 ± 66 
-81 ±65 


-6 ±18 
±3± 18 



Table 2: Results of a study of fragmentation effects on the M w measurement (see text). 
All results are given for the fitted mass analyses on the samples studied with respect to 
the standard Delphi JETSET sample, unless otherwise stated. The statistical error on 
the observed difference is also given. 



T w Fragmentation study (MeV/c 2 ) 


Event Type 


ez7 c qq' 


/^ M qq' 


rV T qq' 


qq'qq' 


Aqcd - 2cr 
Aqcd ± 2(7 

O q - 2(7 

a q ± 2<7 


-5 ±30 
±11 ±30 
±16 ±30 
±28 ± 30 


-46 ± 35 
-47 ±35 
-76 ± 35 
-78 ± 35 


-8 ±59 
±57 ±59 
-11 ±59 
-30 ± 59 


-4± 12 
±15 ± 12 
-11 ± 12 
-2 ± 12 


Aleph HERWIG - AlephjetSET 
Aleph JETSET 


-148 ± 113 
±264 ± 108 


-47 ± 95 
-30 ± 95 


-155 ± 158 
±176 ± 157 


-68 ± 38 
±4 ±43 



Table 3: Results of a study of fragmentation effects on the T-w measurement (see text). 
All results are given for the fitted mass analyses on the samples studied with respect to 
the standard Delphi JETSET sample, unless otherwise stated. The statistical error on 
the observed difference is also given. 



1G 



M w ISR study (MeV/c 2 ) 


Model 


ez/ e qq' /ii^qq' rz/ r qq' 


qq'qq' 


KORALW re-weighted - QEDPS 


+2 ±7 0±5 -12 ±8 


-16 ±3 



Table 4: A comparison of two ISR treatments for the M w measurement, the KORALW form 
and the QEDPS treatment. The statistical error on the result is also provided. 



T w ISR study (MeV/c 2 ) 


Model 


ez/ c qq' /iz^qq' 


rz/ r qq' 


qq'qq' 


KORALW re-weighted - QEDPS 


+13 ± 13 +14 ± 10 


+4+16 


-16 + 5 



Table 5: A comparison of two ISR treatments for the Tyy measurement, the KORALW form 
and the QEDPS treatment. The statistical error on the result is also provided. 



Bose-Einstein Correlations Study (MeV/c 2 ) 


Model 


M w shift 


T w shift 


KKM model (183GeV) 


-10 ± 10 




ST model BE (183GeV) inside W's - none 

ST model BE (183GeV) between W's - inside W's 


+0 ± 10 
+3 ± 11 




LUBOEI BE 32 (189GeV) inside W's - none 
LUBOEI BE m (189GeV) between W's - inside W's 


+18 ± 5 
-32 ± 4 


+ 13 ± 11 
+26 ± 8 



Table 6: Results of the Delphi studies on Bose-Einstein Correlations, see text for details. 
The BE32 samples were produced with A = 1.35, r = 0.6 fm 



Colour Reconnection Study (MeV/c 2 ) 


Model 


M w shift 


r\y shift 


ARIADNE CR 2 (183GeV) 
ARIADNE CR 3 (183GeV) 


+28 ± 6 
+55 ± 6 




JETSET SKI 

JETSET SKI improved sampling 
JETSET SK2 


46 + 2 
44 + 2 
-2 + 5 


54 + 3 
32 + 3 
37 ± 10 



Table 7: Results of the Delphi studies on Colour reconnection effects, see text for details. 
For the SKI model results are given for 30% of reconnected events. 
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M w Systematic Errors (MeV/c 2 ) 



Soiiffp^ of ^v^tpmpfip PT*T*OT* 

kJ w LiX LCO Ul O V ijliClllCHjlU OA 1U1 


P7/ nn' 




TU nn' 


PlJ nC\C\ 


nn nn 
44 44 


Statistical error on calibration 


18 


15 


23 


10 


7 


ijepton collections 


29 


11 




1 n 




Jet corrections 


39 


27 


48 


35 


18 


Aspect Ratio 


2 


2 


2 


2 


4 


Background 


10 


3 


4 


3 


5 


Fragmentation 


20 


20 


20 


20 


12 


I.S.R. 


16 


16 


16 


16 


16 


LEP energy 


17 


17 


17 


17 


17 


Colour reconnection 
Bose Einstein correlations 






46 
32 



Table 8: Contributions to the systematic error on the W mass measurement. 



Tw Systematic Errors (MeV/c 2 ) 



Sources of systematic error 


eT7 c qq' 


/iT7 M qq' 


ru T qq' 


£77 £ qq' 


qq'qq' 


Statistical error on calibration 


45 


37 


56 


26 


17 


Lepton corrections 


41 


46 




28 




Jet corrections 


82 


43 


102 


63 


26 


Background 


29 


8 


82 


19 


40 


Fragmentation 


42 


42 


42 


42 


24 


I.S.R. 


16 


16 


16 


16 


16 


Colour reconnection 
Bose Einstein correlations 






54 
26 



Table 9: Contributions to the systematic error on the W width measurement. 



189 GeV M w results (GeV/c 2 ) 


Channel 


M w stat. syst. LEP FSI 


ez7 e qq' 

/iz7 M qq' 

r^qq' 


80.478 ±0.291 ±0.059 ±0.017 
80.195 ±0.213 ±0.042 ±0.017 
80.114 ±0.319 ±0.059 ±0.017 


£z/£qq' 


80.253 ±0.151 ±0.046 ±0.017 


qq'qq' 


80.466 ±0.106 ±0.028 ±0.017 ±0.056 



Table 10: 189 GeV Mw results. The error is divided into its statistical component, 
indicated stat, the main systematic component syst and the systematic from the beam 
energy uncertainty LEP. In the qq'qq' channel an error from final state interaction 
effects FSI is also included. The £^qq' results represents the combination of the results 
obtained in the three semi-leptonic channels. 
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189 GeV T w results (GeV/c 2 ) 


Channel 


Tw stat. syst. FSI 


ez7 e qq' 
rz7 T qq' 


4.358 ±0.956 ±0.115 
2.353 ±0.552 ±0.086 
2.799 ±0.927 ±0.149 




2.842 ±0.425 ±0.088 


qq'qq' 


2.025 ±0.220 ±0.058 ±0.060 



Table 11: 189 GeV T\y results. The error is divided into its statistical component, 
indicated stat, the main systematic component syst and the systematic from the beam 
energy uncertainty LEP. In the qq'qq' channel an error from final state interaction 
effects FSI is also included. The iu^qq' results represents the combination of the results 
obtained in the three semi-leptonic channels. 




-0.2 0.2 0.4 0.6 0.8 1 1.2 

neual netwok output br qq xv candidates 



Figure 1: The output distribution of the event selection neural network for rz7 r qq' candi- 
dates. The shaded areas indicate the contribution from the various simulated states, the 
data are shown as points with statistical error bars. The value of the selection cut applied 
is indicated by the dashed line. Note that the order of the simulation contributions in 
the figure follows that in the key, the rV T qq' signal is the first distribution shown and 
background from two-photon diagrams the last. 
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(a) 



(b) 



Figure 2: An example of the reconstructed probability density 

function J2i w i,e • Pi,e(m+, m -) (see text) for the invariant masses in a simulated 4-jet 
qq'qq' event without (a) and with (b) the hypothesis of collinear ISR. The first 3 sigma 
contours are shown. The normalization of the different solutions prevents the high mass 
contours from reaching the 1 sigma probability level, while the small difference in the low 
mass solutions originates from the jet charge information. The generated masses of the 
two W bosons in the event are marked with a cross. 
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W mass (GeV/c") 
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(c) 



m 70 
u 
> 

- 60 
ri 

3 50 
e 

40 
30 
20 
10 



> 

o 

O 



WW qqeu 

• data 

qqei/ simulation 
(m w =80.35) 

| qqri/ simulotion 
I Z(y) + ZII simulation 




50 55 60 65 70 75 80 85 90 95 100 



W mass (GeV/c ) 



(b) 
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Figure 3: The distribution of the reconstructed W masses from a kinematic fit with five 
constraints imposed in the (a) qq'qq' , (b) e/7 c qq', (c) /iF^qq' and (d) rF T qq' analysis 
channels. In the qq'qq' channel, only the jet pairing with the highest probability is 
included in this figure. 
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reconnected fraction (%) reconnected fraction (%) 



Figure 4: Observed shift in the fitted mass (left) and width (right) as a function of the 
fraction of reconnected events, using the SKI implementations as described in the text. 



